Overcoming the coherence time barrier in quantum machine learning on temporal data

The practical implementation of many quantum algorithms known today is limited by the coherence time of the executing quantum hardware and quantum sampling noise. Here we present a machine learning algorithm, NISQRC, for qubit-based quantum systems that enables inference on temporal data over durations unconstrained by decoherence. NISQRC leverages mid-circuit measurements and deterministic reset operations to reduce circuit executions, while still maintaining an appropriate length persistent temporal memory in the quantum system, confirmed through the proposed Volterra Series analysis. This enables NISQRC to overcome not only limitations imposed by finite coherence, but also information scrambling in monitored circuits and sampling noise, problems that persist even in hypothetical fault-tolerant quantum computers that have yet to be realized. To validate our approach, we consider the channel equalization task to recover test signal symbols that are subject to a distorting channel. Through simulations and experiments on a 7-qubit quantum processor we demonstrate that NISQRC can recover arbitrarily long test signals, not limited by coherence time.

The development of machine learning algorithms that can handle data with temporal or sequential dependencies, such as recurrent neural networks [1] and transformers [2], has revolutionized fields like natural language processing [3].Realtime processing of streaming data, also known as online inference, is essential for applications such as edge computing, control [4], and forecasting [5].The use of physical systems whose evolution naturally entails temporal correlations appears at first sight to be ideally suited for such applications.An emerging approach to learning, referred to as physical neural networks (PNNs) [6][7][8][9], employs a wide variety of physical systems to compute a trainable transformation on an input signal.A branch of PNNs that has proven well suited to online data processing is physical reservoir computing [10], distinguished by its trainable component being only a linear projector acting on the observable state of the physical system [11].This approach has the enormous benefit of fast convex optimization through singular value decomposition routines and has already enabled temporal learning on various hardware platforms [4,[12][13][14][15].Among many physical systems considered for PNNs, quantum systems are believed to offer an enormous potential for more scalable, resource-efficient, and faster machine learning [16][17][18][19][20][21][22][23], due to their evolution taking place in the Hilbert space that scales exponentially with the number of nodes [24][25][26][27][28][29][30].However, quantum machine learning (QML) on presentday noisy intermediate-scale quantum (NISQ) hardware has so far been restricted to training and inference on lowdimensional static data due to several difficulties.A fundamental restriction is Quantum Sampling Noise (QSN) -the unavoidable uncertainty arising from the finite sampling of a quantum system -which limits the accuracy of both QML training and inference [9,31,32] even on a fault-tolerant hardware.Additionally, the optimization landscape for training quantum systems often features "barren plateaus" [33,34], which are regions where optimization becomes exponentially difficult.These plateaus, especially in the presence of QSN, present a significant challenge to implementing QML at scales relevant to practical applications.Two further concerns arise when considering inference on long data streams, which call into question whether quantum systems can even in principle be employed for online learning on streaming data.Firstly, without quantum error correction, the operation fidelities and finite coherence times of constituent quantum nodes places a limit on the size of data on which inference can be performed [35,36], which would appear to rule out inference on long data streams.Secondly, the nature of measurement on quantum systems imposes a fundamental constraint on continuous information extraction over long times.Backaction due to repeated measurements on quantum systems necessitated by inference on streaming data is expected to lead to rapid distribution of information between different parts of the system, a phenomenon known as information scrambling and thermalization [37,38], making it extremely difficult to track or retrieve the information correlations in the input data.This constraint persists even in an ideal system with perfect coherence, such as one that may be realized by a fault-tolerant quantum computer.It is not known precisely what conditions must be satisfied to avoid information scrambling.For classical dynamical systems, a strict condition known as the fading memory property [39,40] is required for a physical system to retain a persistent temporal memory that does not degrade on indefinitely long data streams.This imposes restrictions on the design of a classical reservoir and in particular how input data is encoded.Here, a mathematical framework known as Volterra Series theory [41] provides the basis for analyzing the memory properties of a classical dynamical system.Such a general theory for quantum systems has remained elusive so far.
Here we present a Volterra theory for quantum systems that accounts for measurement backaction, necessary for analyzing the conditions required to endow a quantum system with a persistent temporal memory on streaming data.Based on this Quantum Volterra Theory we propose an algorithm, NISQ Reservoir Computing (NISQRC), that leverages recent technical advances in mid-circuit measurements to process signals of arbitrary duration, not limited by the coherence time of constituent physical qubits (see Fig. 1).The property that enables inference on an indefinitely-long input signal -the ability to avoid measurement-induced thermalization at long times under repeated measurements due to a deterministic reset protocol -is intrinsic to the algorithm: it survives even in the presence of QSN, and does not require operating in a precisely-defined parameter subspace -and is thus unencumbered by barren plateaus.
Here, we demonstrate the practical viability of NISQRC through application to a task of technological relevance for communication systems, namely, the equalization of a wireless communication channel.Channel equalization aims to reconstruct a message streamed through a noisy, non-linear and distorting communication channel and has been employed in benchmarking reservoir computing architectures [11,14] as well as other machine learning algorithms [42,43].This task poses a challenge for parametric circuit learning-based algorithms [19] because the number of symbols in the message, N ts , to recover in the inference stage directly determines the length of the encoding circuit, which in turn is limited by the coherence time of the system.A more critical issue is that the recovery has to be done online, as the message is streamed, which structurally is not suitable for static encoding schemes.We demonstrate through numerical simulation (Results' subsection "Practical machine learning using temporal data") and experiments on a 7-qubit quantum processor (Results' subsection "Experimental results on quantum system") that NISQRC enables quantum systems to process signals of arbitrary duration.Most significantly, this ability to continuously extract useful information from a single quantum circuit is not limited by coherence time.Instead, the quantum system's coherence influences the resulting memory timescale; we show that by balancing the length of individual input encoding steps with the rate of information extraction through mid-circuit measurements, it is possible to endow the circuit with a memory that is appropriate for the ML task at hand.Even in the limit of infinite coherence, the temporal memory is still limited by this fundamental trade-off.Reliable inference on a time-dependent signal of duration T run = 117µs is demonstrated on a 7-qubit quantum processor with qubit lifetimes in the range 63µs -164µs and T 2 = 9µs -231µs.In our experiments longer durations are restricted by limitations on mid-circuit buffer clearance.To leave no doubt that a persistent memory can be generated, we first compare the experimental results to numerical simulations with the same parameters, showing excellent agreement.Building on the accuracy of numerical simulations in the presence of finite coherence  1.Schematic representation of NISQRC architecture for machine learning on temporal data using a convex optimization algorithm on finitely-sampled partial measurements.For concreteness, the architecture is shown for a quantum circuit with projective computational basis readout; both the underlying quantum system and the measurement scheme can be much more general.Temporal input data is encoded into the evolution of the reservoir at every time-step n via a quantum channel U(un); a non-trivial I/O map is enabled via partial readout and subsequent reset of a readout subsystem, while a memory subsystem retains memory of past inputs.Temporal quantum reservoir computing (QRC) output x(n) are obtained via measurements (more precisely, stochastic unbiased estimators X(n) of expected features are constructed from S repetitions of the experiment, see Method III A), and a learned linear combination is used to approximate the target functional y(n) of un.The overall execution time of the circuit is O(N S), where N is the length of input temporal sequence.and our noise model, we explicitly demonstrate successful inference on a 5000 symbol signal: the resulting circuit duration is 500 times that of the individual qubit lifetimes.
Here, we also develop a method to efficiently sample from deep circuits under partial measurements.Simulating individual quantum trajectories for circuits with repeated measurements requires the traversal of ever-branching paths conditioned on the measurement results, which becomes rapidly unfeasible for deep circuits.Our numerical method (see Methods' subsection "The quantum Volterra theory and analysis of NISQRC") allows us to sample from repeated partial measurements on circuits of arbitrary depth.We use our scheme to numerically explore other seemingly reasonable encoding methods adopted in previous studies, showing that these can lead to a sharp decline in performance when the effect of measurement is properly accounted for.Drawing upon the Quantum Volterra Theory, we unveil the underlying cause: the absence of a persistent memory mechanism.

A. Time-series processing in quantum Systems
The general aim of computation on temporal data is expressed most naturally in terms of functionals of a time- A functional F : u → y maps a bounded function u to another arbitrary bounded function y, where y Without loss of generality these functions can be normalized; we choose u n ∈ [−1, 1] and y n ∈ [−1, 1].Within the reservoir computing paradigm [44], this processing is achieved by extracting outputs x(n), where n is a temporal index, from a physical system evolving under said time-dependent stimulus u n ≡ u(n).Learning then entails finding a set of optimal timeindependent weights w to best approximate a desired F with a linear projector y n ≡ y(n) = w • x(n).If the physical system is sufficiently complex, its temporal response x(n) to a time-dependent stimulus u is universal in that it can be used to approximate a large set of functionals F[u] with an error scaling inversely in system size and using only this simple linear output layer [27,28,45].
To analyze the utility of this learning framework, it proves useful to quantify the space of functionals F[u] that are accessible.For classical non-linear systems, a firmly-established means of doing so is a Volterra series representation of the input-output (I/O) map [39]: (1) where the Volterra kernels h k over the the temporal domain (n 1 , • • • , n k ) quantifies the notion of memory of a particular physical system, with the kernel order k being the corresponding degree of nonlinearity of the map.Most importantly, the Volterra series representation describes a time-invariant I/O map, as well as the property of fading memory, which roughly translates to the property that the reservoir forgets initial conditions and thus depends more strongly on more recent inputs (for instance, for multi-stable dynamical systems, a global representation such as Eq. ( 1) may not exist.However a local representation around each steady state can be shown to exist with a finite convergence radius).The realization of such a time-invariant map is essential for a physical system to be reliably employed for inference on an input signal of arbitrary length, and thus for online time series processing.
In classical physical systems, the existence of a unique information steady state and the resulting fading memory property is determined only by the input encoding dynamics -the map from input series to system state.More explicitly, the information extraction step (sometimes referred to as the "output layer") on a classical system is considered to be a pas-sive action, so that the state can always be observed at the precision required.However, for physical systems operating in the quantum regime, the role of quantum measurement is fundamental: in addition to the inherent uncertainty in quantum measurements as dictated by the Heisenberg uncertainty principle, the conditional dependence of the statistical system state on prior measurement outcomes -referred to as backaction -strongly determines the information that can be extracted.Recent work in circuit-based quantum computation has shown that the qualitative features of the statistical steady state of monitored circuits strongly depends on the rate of measurement [46,47].In particular, generic quantum systems that alternate dynamics and measurement (input encoding and output in the present context) are known to give rise to deep thermalization of the memory subsystem [48,49], resulting in an approximate Haar-random state with vanishing temporal memory.The absence of a comprehensive framework in QML for analyzing and implementing an encoding-decoding system with finite temporal memory, along with characterization tools for the accessible set of input-output functionals, has hindered both a systematic study and the practical application of online learning methods.
Here, we develop both a general temporal learning framework suitable for qubit-based quantum processors and the associated methods of analysis based on an appropriate generalization of the Volterra Series analysis to monitored quantum systems, the Quantum Volterra Theory (QVT).Our approach incorporates the effects of backaction that results from quantum measurements in the process of information extraction.
We begin by providing a fundamental description of both the information input and output processes that enable general time series processing with quantum systems, before specializing to the NISQRC algorithm.The 'input' component of the map is given by a pipeline (encoding) that injects temporal data {u n } into a quantum system through a general parameterized quantum channel U(u n )ρ .This channel could describe for instance continuous Linblad evolution for a duration τ , namely e τ L(un) ρ, as in Results' subsection "Practical machine learning using temporal data", or a discrete set of gates as in Results' subsection "Experimental results on quantum system"; U(u n ) is generally applied to all qubits, and we assume only that they are not explicitly monitored for its duration.
To enable persistent memory in the presence of quantum measurement, we separate the L-qubit system into M memory qubits and R readout qubits (L = M + R), and denote their respective Hilbert spaces with superscript M and R.After evolution under any input u n , only the R readout qubits are (simultaneously) measured; this separation therefore allows for the concept of partial measurements of the full quantum system, which proves critical to the success of NISQRC.The measurement scheme itself can be very general, characterized by a positive operator-valued measure (POVM) satisfying Êj ⪰ 0 and j Êj = Î⊗R .Here we will con-sider a practically implementable measurement in the readout qubit computational basis, described by Êj = |b j ⟩⟨b j |: each bit-string b j is the R-bit binary representation of integer j ∈ {0, 1, • • • , 2 R − 1} denoting the bit-wise state of the measured qubits.
As elucidated by the QVT analysis of Results' subsection "Quantum Volterra Theory", a purification mechanism must necessarily accompany readout to prevent thermalization and furnish our quantum architecture with persistent fading memory.This is accomplished by following each projective measurement operation with a deterministic reset to the ground state |0⟩.The resulting measure-reset operation we employ throughout this paper is formally described by the POVM operators Êj = K † j Kj in Eq. ( 2), with non-hermitian Kraus operators Kj = |b 0 ⟩⟨b j |.In each measure-reset step, only the readout qubits are measured in the computational basis and then reset to the ground state, irrespective of the measurement outcome.
NISQRC is distinguished by the iterative encode-measurereset scheme depicted in Fig. 1.Explicitly, for a given input sequence u with length N , we initialize the system in the state ρM 0 ⊗ |0⟩⟨0| ⊗R .For each element of the input sequence u n , an encoding step is comprised of unmonitored evolution of all qubits via U(u n ) followed by a measurereset operation O R .The measurement outcome in this single shot is a random bitstring b (s) (n), and the resulting state is ρM,cond n

⊗ |0⟩⟨0|
⊗R : the memory qubits are in a state conditioned on the measurement outcome, and the readout qubits are reset.The subsequent input is then encoded in this state, i.e.U(u n+1 ) ρM,cond n

⊗|0⟩⟨0|
⊗R , and the process is iterated as long as there is data in the pipeline.This structure elucidates the naming of the unmeasured memory qubits: these are the only qubits that retain memory of past inputs.
The above description yields a set of N measurement outcomes {b (s) (n)} observed in a single shot s of the quantum circuit.In order to obtain statistics and therefore to output features as expected values of observables Mj , we perform S repetitions of this circuit for a given u sequence: the total execution time is N S, linear with respect to shots S and input length N .The resulting readout features are formally defined as the probability which are estimated by the empirical mean (see Methods' subsection "Generating features via conditional evolution and measurement" for more details of NISQRC algorithm).We show in Supplementary Note 2 that at time step n, x j (n) can be computed efficiently through where ρMR n is the effective full L-qubit system state at time step n prior to measurement.
The output y n ≡ y(n) = w • x(n) ∈ R is obtained from the measurement results in each step, defining the functional I/O map which we characterize next (see details in Methods' subsection "Generating features via conditional evolution and measurement" and "The quantum Volterra theory and analysis of NISQRC").This complete architecture, from the quantum circuit generating measurement outcomes for a given input, to the construction of weighted output features, is depicted schematically in Fig. 1.We note that reset operations have been used implicitly in prior work on quantum reservoir algorithms, where the successive inputs are encoded in the state of an 'input' qubit [26,50].However the critical role of the reset operation in endowing a quantum reservoir with a persistent memory, discussed in the next section, has so far not been highlighted.While for null inputs (i.e.u n = 0 for all n) such quantum systems are guaranteed to have a unique statistical steady state, the existence of a nontrivial memory and kernel structure is much more involved.Through QVT (see Methods' subsection "The quantum Volterra theory and analysis of NISQRC"), we show that these requirements place strong constraints on the encoding and measurement steps viz. the choice of (U, Mj ).This then enables us to propose an algorithm for online learning that provably provides a controllable and time-invariant temporal memory (which will be referred to as persistent memory) -enabling inference on arbitrarily long input sequences even on NISQ hardware without any error-mitigation or correction.

B. Quantum Volterra theory
In NISQRC the purpose of the partial reset operation is to endow the system with asymptotic time-invariance, a finite persistent memory and a nontrivial Volterra Series expansion for the system state (see Methods' subsection "The quantum Volterra theory and analysis of NISQRC" and Supplementary Note 3): where all Volterra kernels ĥk are quantum operators.The classical kernels in Eq. ( 1) describing the measured features can be extracted through h (j) k = Tr( Mj ĥk ).We refer to this analysis as the Quantum Volterra Theory (QVT).Through analytical arguments based on the QVT, we show that omitting the partial reset operation renders all Volterra kernels trivial -a finding corroborated by our experimental results in Results' subsection "Experimental results on quantum system".
QVT also provides a way to characterize the important memory time-scales of the I/O map generated by the NISQRC algorithm through a given encoding, which we use in Results' subsection "Practical machine learning using temporal data" to aid encoding design for a specific ML task on an experimental system.In what follows, we show that inference on an indefinitely long input sequence can be done even in the presence of dissipation and decoherence.FIG. 2. Quantum Volterra Theory (QVT) analysis for (M + R)qubit reservoir.(a) First and second order Volterra kernels in a (2 + 1)-qubit QRC, which vanish at large n1 and n2 due to finite memory nM.(b) Fixed-point of memory subsystem ρM FP with reset (top) and without reset (bottom), starting from an arbitrary initial state (center).Without reset, the fixed point is always the trivial fully-mixed state and Volterra kernels vanish.Top panel shows the distribution of the 4 M = 256 eigenvalues of P0 in a (4 + 2)-qubit QRC, where red dots correspond to the static unit eigenvalue λ1 = 1.The remaining eigenvalues λ α≥2 (blue) evolve with evolution time τ , leading to a variable memory time.Bottom panel shows the resulting memory time nM as a function of the evolution duration τ .(c) Memory time nM as a function of qubit lifetimes T1 = γ −1 , in terms of the evolution duration τ in a (4 + 2)-qubit QRC.Provided T1 ≫ τ , nM → n 0 M , so that the QRC memory is mostly dominated by its lossless dynamical map, and not by T1 in this regime.
Consider an input-encoding U(u n )ρ = e τ L(un) ρ where representing evolution under a parameterized Hamiltonian Ĥ(u n ) for a duration τ in the presence of dissipation D T .
For concreteness, we take ] describing decoherence processes and study here a specific Ising Hamiltonian encoding Ĥ(u) = Ĥ0 + u • Ĥ1 inspired by quantum annealing and simulation architectures (other ansätze can likewise be considered), The coupling strength J i,i ′ , transverse x-field strength η x i and longitudinal z-drive strength η z i are randomly chosen, but then fixed for all inputs {u n } (see Supplementary Note 1 for more details).The encoding channel is applied for duration τ , and each qubit has a finite lifetime T 1 = γ −1 .We will specify the number of memory and readout qubits of a given QRC with the notation (M + R).
In Fig. 2(a) we plot the first two Volterra kernels h 1 and h 2 (cf.Eq. ( 1)) for a random (2 + 1)-qubit QRC using the above encoding and the reset scheme.The expression for these kernels have been derived from the QVT; their numerical construction is discussed in Methods, also see Supplementary Equations 43-45.Importantly, we find all kernels have an essential dependence on the statistical steady state or fixedpoint in the absence of any input: ρM FP = lim n→∞ ρM n un=0 .Here ρM n un=0 = P n 0 ρM 0 is obtained by n applications of the null-input single-step quantum channel P 0 , defined in Methods' subsection "The quantum Volterra theory and analysis of NISQRC".The properties of quantum Volterra kernels, including their characteristic decay time, can be related to the spectrum of P 0 , defined by P 0 ρM α = λ α ρM α .Here ρM α are eigenvectors that exist in the 4 M -dimensional space of memory subsystem states.The eigenvalues satisfy The second largest eigenvalue λ 2 determines the time over which memory of an initial state persists as this fixed point is approached, and is used to identify a memory time n M = −1/ ln |λ 2 |.Note that this quantity is dimensionless and can be converted to actual passage of time through multiplication by τ , while n M itself non-trivially depends on τ (see Fig. 2(b)).The memory time describes an effective 'envelope' for a system's Volterra kernels; additional nontrivial structure is also required for QRC to produce meaningful functionals of past inputs.With the spectral problem at hand, we next analyze the information-theoretical benefit of the reset operation.Firstly, the absence of the unconditional reset operation produces a unital P 0 ("unital" refers to an operator that maps the identity matrix to itself).with resulting ρM FP = I ⊗M /2 M .This fully-mixed state is inexorably approached after n M steps under any input sequence and retains no information on past inputs: all Volterra kernels therefore vanish, despite a generally-finite n M .Such algorithms (e.g.Refs.[27,30]) are only capable of processing input sequences of length n M and would not retain a persistent memory necessary for inference on longer sequences of inputs.Hence such encodings would be unsuitable for online learning on streaming data.The possibility of inference through the transients have been observed and utilized before (see e.g.Ref. [12,51,52]) in the context of classical reservoir computing.However, the simple yet essential inclusion of the purifying reset operation avoids unitality -more generally, a common fixed point for all u-encoding channels -which we find is the key to enabling nontrivial Volterra kernels and consequent online QRC processing (see Methods' subsection "The quantum Volterra theory and analysis of NISQRC" and also Ref. [53]).Once such an I/O map is realized, λ α and the consequent memory properties can be meaningfully controlled by the QRC encoding parameters.As shown in Fig. 2(b) the characteristic decay time set by n M , for instance, decreases across several orders-of-magnitude with increasing τ .
The partial measurement and reset protocol also resolves the unfavorable quadratic runtime scaling of prior approaches.
A wide range of proposals and implementations of QRC [27,29,54] consider the read out of all constituent qubits at every output step, terminating the computation.Not only does this preclude inference on streaming data, it requires the entire input sequence to be re-encoded to proceed one step further in the computation, leading to an O(N 2 S) running time.As shown in schematic Fig. 1, incorporating partial measurement with reset in NISQRC does not require such a re-encoding; the entire input sequence can be processed in any given measurement shot S, enabling online processing with an O(N S) runtime, while maintaining a controllable memory timescale.We note that an alternative scheme to remedy this issue has been suggested in Ref. [50], which relies on information extraction through continuous weak measurement.
Next we show that the nontrivial nature of Volterra kernels realized by the NISQRC algorithm is preserved under the inclusion of dissipation.For example, we explore the effect of finite qubit T 1 on n M in Fig. 2(c).If T 1 /τ > n 0 M , where n 0 M is the memory time of the lossless map, then n M → n 0 M and is essentially independent of T 1 , determined instead by the unitary and measurement-induced dynamics.Therefore the design of the encoding algorithm has to be guided by matching the memory time of the reservoir to the longest correlation time in the input data.Additional design criteria are discussed in Section Discussion.As a result, lossy QRCs can still be deployed for online processing, with a total run time T run that is unconstrained by (and can therefore far exceed) T 1 .We will demonstrate this via simulations in Results' subsection "Practical machine learning using temporal data" with T run ≫ T 1 , and via experiments in Results' subsection "Experimental results on quantum system" for T run ≃ T 1 ; in the latter T run is limited only by memory buffer constraints on the classical backend.

C. Practical machine learning using temporal data
Thus far, we have assumed outputs to be expected features x j (n), which in principle assumes an infinite number of measurements.In any practical implementation, one must instead estimate these features with S shots or repetitions of the algorithm for a given input u.The resulting QSN constrains the learning performance achievable in experiments on quantum processors in a way that can be fully characterized [9], and is therefore also included in numerical simulations which we present next.
To demonstrate the utility of the NISQRC framework, we consider a practical application of machine learning on time-dependent classical data: the channel equalization (CE) task [11,14].Suppose one wishes to transmit a message m(n) of length N , which here takes discrete values from {−3, −1, 1, 3}, through an unknown noisy channel to a receiver.This medium generally distorts the signal, so the received version u(n) is different from the intended m(n).Channel equalization seeks to reconstruct the original message m(n) from the corrupted signal u(n) as accurately as Without reset (red) or using 4 ancilla qubit ansatz with quantum non-demolition (QND) readout (proposed in Ref. [30], green), the algorithms both fail, approaching the random guessing error rate and showing that both architectures suffer from the thermalization problem.Performance is only slightly reduced from the dissipationfree case (blue) when strong decay T1 = 10τ is included (purple).All error rates in (c) are averaged over 8 different test messages.possible, and is of fundamental importance in communication systems.Specifically, we assume the message is corrupted by nonlinear receiver saturation, inter-symbol interference (a linear kernel), and additive white noise [11,14] (additional details in Supplementary Note 6).As shown in Fig. 3(a), even if one has access to the exact inverse of the resulting nonlinear filter, the signal-to-noise (SNR) of the additive noise bounds the minimum achievable error rate.We also show the error rates of simply rounding u(n) to the nearest m, and a direct logistic regression on u(n) (i.e. a single-layer perceptron with a softmax activation -see Supplementary Note 6).for comparison.Both these approaches are linear and memory-less and therefore perform poorly on the non-trivial nonlinear filter we consider, although logistic regression outperforms rounding (≈30%) by inverting the linear portion of the distortion.
We now perform the CE task using the NISQRC algorithm on a simulated (2 + 4)-qubit reservoir under the ansatz of Eq. ( 7), as could be realized in quantum annealing hardware (see Fig. 3).We will later demonstrate the same task in experiments with a completely different quantum system and encoding ansatz, implemented on a superconducting quantum processor (see Fig. 4).The ability to efficiently compute the Volterra kernels for this quantum system immediately provides guidance regarding parameter choices.In particular, we choose random parameter distributions such that the average (across the circuit) J i,i ′ τ , η x i τ and η z i τ provides a memory time n M ≈ O(10 1 ), on the order of the length of the distorting linear kernel h(n), which is 8.These QRCs have K = 2 4 = 16 readout features {x j (n)} j∈[K] whose corresponding time-independent output weights w are learned by minimizing cross-entropy loss on 100 training messages of length N = 100 (see Supplementary Note 6 for additional details).The resulting NISQRC performance on test messages is studied in Fig. 3(a), where we compare two distinct coupling maps shown in Fig. 3(b).In the highly-connected (lower) system the performance approaches the theoretical bound for S → ∞; finite sampling (here, S = 10 5 is in the range typically used in experiments) increases the error rate as expected, but the increase in error rate in numerical simulations is observed to depend on the encoding (not reported here).In all cases, NISQRC significantly outperforms direct logistic regression due to its ability to reliably implement nonlinear memory kernels and therefore approximate the distorting channel inverse.
We note that the split system (upper) performs significantly worse even without sampling noise: this is because the quantum system lives in a smaller effective Hilbert space -the product of two disconnected three-qubit systems -and is far less expressive as a result.Although in both cases the number of measured features is the same, those from the connected system span a richer and independent space of functionals.This functional independence can be quantified by the Jacobian rank R J , which is the number of independent u-gradients that can be represented by a given encoding (Supplementary Note 5); an increased connectivity and complexity of statedescription generally manifests as an increase in the Jacobian rank and consequent improved CE task performance.This observation can be viewed as a generalization of the findings in time-independent computation [9] to tasks over temporallyvarying data, and also agrees with related recent theoretical work [29].
Most importantly, we demonstrate in Fig. 3(c) that the NISQRC algorithm enables the use of a quantum reservoir for online learning.In all cases studied here, N = 100 is used for training and the length of the SNR = 20 dB test messages N ts is varied.As suggested by the QVT, the performance is unaffected by N ts even if it greatly exceeds the lifetime of individual qubits: N ts = T run /τ ≫ T 1 /τ = 10, and NISQRC can therefore be used to perform inference on an indefinite-length signal with noisy quantum hardware.As seen in the same figure, while dissipation imposes only a small constant performance penalty, the reset operation is critical: if removed, the error rate increases to that of random guessing, as the Volterra kernels vanish and the I/O map becomes trivial.
We finally note that an arbitrarily-inserted reset operation may not be sufficient to create a non-zero persistent memory.For instance, an analysis based on the QVT shows that despite its use in a recently studied reservoir algorithm [30] (based on a quantum non-demolition measurement proposal in Ref. [27]), the reset operation can not avoid a zero persistent memory, effectively resulting in an amnesiac reservoir.In this scheme, the quantum circuit is coupled to ancilla qubits by using transversal CNOT gates.Upon closer examination it is found that while the projective measurement of ancilla qubits leads to read out of system qubits and their collapse to the ancilla state via backaction, subsequent reset of the ancillas does not reset the system qubits.This scheme therefore suffers from the same thermalization problem as any no-reset NISQRC does, and hence has zero persistent memory.We verify this analysis in Fig. 3(c) by implementing the CE task with a four-ancilla-qubit circuit.The error rates are found to be very close to the no-reset-NISQRC one, whose I/O map we have shown before to be trivial (see also Fig. 3(c)).

D. Experimental results on quantum systems
We now demonstrate NISQRC in action by performing the SNR = 20 dB CE task on an IBM Quantum superconducting processor.To highlight the generality of our NISQRC approach, we now consider a circuit-based parametric encoding scheme inspired by a Trotterization of Eq. ( 7), suitable for gate-based quantum computers.In particular, we use a L = 7 qubit linear subgraph of the ibm algiers device, with M = 3 memory qubits and R = 4 readout qubits in alternating positions, as depicted in Fig. 4(a).The encoding unitary for each time step n is also shown: n T , where R x,z are composite Pauli-rotations applied qubit-wise, and W(J) defines composite R zz gates between neighbouring qubits, all repeated n T = 3 times (for parameters θ x,z,I , J and further details see Methods' subsection "IBMQ implementation").
Realizing the NISQRC framework with the circuit ansatz depicted in Fig. 4(a) requires the state-of-the-art implementation of mid-circuit measurements and qubit reset, which has recently become possible on IBM Quantum hardware [55].We plot the testing error using the indicated linear chain of the ibm algiers device as a function of the number of shots S in solid blue Fig. 4(b), alongside simulations of both the ideal unitary circuit and with qubit losses in open circles.We clearly observe that performance is influenced by the number of shots available, and hence by QSN.In particular, for a sufficiently large S, the device outperforms the same logistic regression method considered previously.For the circuit runs, the average qubit coherence times over 7 qubits are T av 1 = 124 µs, T av 2 = 91µs (see Supplementary Note 9 for the ranges of all parameters, which varies over the time of runs as well), while the total circuit run time for a single message is T run ≈ 117 µs.Even though T run ≃ T av 1 , the CE task performance using NISQRC on ibm algiers is essentially independent of qubit lifetimes.This is emphatically demonstrated by the excellent agreement between the experimental results and simulations assuming infinite coherence-time qubits.In fact, finite qubit decay consistent with ibm algiers leaves simulation results practically unchanged (as plotted in dashed blue); we find that T 1 times would have to be over an order of magnitude shorter to begin to detrimentally impact NISQRC performance on this device (see Supplementary Note 7).We further find that artificially increasing T run beyond T 1 by introducing controlled delays in each layer also leaves performance unchanged (see Supplementary Note 8).
Using the same device we are able to analyze several important aspects of the NISQRC algorithm.First, we consider the same CE task with a split chain, where the connection between the qubits labelled '14' and '16' on ibm algiers is severed by removing the R zz gate highlighted in brown in Fig. 4(a).The resulting device performance using these two smaller chains is worse, consistent both with simulations of the same circuit and the analogous split Hamiltonian ansatz studied in Results' subsection "Practical machine learning using temporal data".Next we return to the 7 qubit chain but now remove reset operations in the NISQRC architecture, shaded in red in Fig. 4(a): all other gates and readout operations are unchanged.The device performance now approaches that of random guessing: the absence of the crucial reset operation leads to an amnesiac QRC with no dependence on past or present inputs.This remarkable finding reinforces that reset operations demanded by the NISQRC algorithm are therefore essential to imbue the QRC with memory and enable any nontrivial temporal data processing.
We note that for these experiments, while performance qualitatively agrees well with simulations, some quantitative discrepancies are observed.Our deployment of mid-circuit measurements in their earliest implementation on IBM Quantum were accompanied by some technical constraints; for example, not all shots for a given instance of the CE task could be collected in contiguous repeated device runs, instead sometimes being separated by several hours (due to queuing times as well as classical memory buffer constraints on the number of shots that could be collected in a single experiment).Simply put, this means that the device could suffer non-trivial parameter drifts from one type of device configuration to the next, and even during the course of collecting all shots for a specific configuration.In particular, we find that qubit lifetimes for experiments with the split chain, and the connected chain without reset, were significantly shorter than for the connected chain with reset (see Supplementary Tables 1-3), which could lead to the discrepancy in comparison to simulations, where we assumed a fixed coherence time distribution.Resource constraints similarly restrict us to limited training and testing set sizes, which can also lead to variance in performance.We anticipate such technical constraints to be alleviated as mid-circuit measurement implementations mature on IBM Quantum, enabling even more accurate correspondence with simulations.
We also note that there is room for improvement in CE performance when compared against Hamiltonian ansatz NISQRC of similar scale in Fig. 3.A key difference is the reduced number of connections in the nearest-neighbour linear chain employed on ibm algiers; including effective R zz gates between non-adjacent qubits significantly increases the gatedepth of the encoding step, enhancing sensitivity to circuitfidelity due to increasing runtimes.The circuit ansatz can also be optimized -using knowledge of the Volterra kernels -for better nonlinear processing capabilities demanded by the CE task, in addition to memory capacity determined by n M .Nonetheless, the demonstrated performance and robustness of the NISQRC framework to dissipation already suggests its viability for increasingly complex time-dependent learning tasks using actual quantum hardware.

II. DISCUSSION
A key technical advancement in our work is the formulation of the Quantum Volterra Theory (QVT) to describe the timeinvariant input-output map of a quantum system under temporal inputs and repeated measurements.Insights provided by the QVT enabled us to propose the essential component of the NISQRC algorithm -deterministic post-measurement reset to avoid thermalization due to repeated measurementswhich allows the quantum system to retain persistent memory of temporal inputs even under projective measurements and their associated strong backaction.The resulting algorithm enables inference on a signal that can be arbitrarily long, provided the encoding is designed to endow the reservoir with a memory that matches the longest correlation time in the data.
While we have applied the QVT to qubit-based circuits, our analysis does not make an explicit assumption on the Hilbert space dimension of the quantum system, and allows for completely general measurements through its formulation in terms of POVMs; as a result, it can be applied to other finite-level quantum systems such as qudits [24], and can be extended to continuous-variable quantum systems [28,31].We therefore believe the QVT provides the ideal framework to analyze the memory and computational capacity of temporal information processing schemes using general quantum systems and their associated measurement protocols.We note here that the use of continuous weak measurements, analyzed in Ref. [50], provides an alternative approach to endowing the reservoir with a finite persistent memory and can be analyzed with QVT for its task-specific optimization.
Going beyond the crucial reset component, we have demonstrated that QVT can be invaluable in identifying general design principles for qubit-based systems as reservoirs.For example, while measuring some fraction of qubits is essential for extracting information, measuring all qubits imposes a trivial memory time.We employ M ≃ R in this work, but an optimal separation of memory and readout qubits may depend on specific tasks.A simple rule of thumb is to choose M/R, together with τ and other drive strengths, to match the memory time of the physical system to the longest correlation time in the data.In the channel equalization task studied in Results' subsection "Practical machine learning using temporal data", the data correlation time is fixed by the choice of the distorting channel and we have then chosen M/R to endow the quantum system with a memory time -calculated through the QVT formalism -that matches that time, about 8 steps (recall memory time is measured in number of encoding steps).Especially, the duration of the unmonitored dynamics, τ , has been chosen to be long enough to generate non-linear kernels that match the known order of the non-linearity of the distorting channel, but short enough to avoid limiting the memory time by the shortest T 1 .For the latter requirement, the kind of analysis shown in Fig. 2(c), calculated through QVT, can act as a very helpful guide.We have observed that even when there is a large spread in T 1 , the physical memory time may be longer than the shortest T 1 , presumably through the delocalization of the information on longer-lived memory qubits in the circuit.Ad-ditionally, qubit connectivity, analyzed in Fig. 3(a), can help with the generation of functions that are sufficiently complex to match the functional complexity of the task.
Finally, the most crucial design criterion for any quantum system intended to process streaming data is that the map P 0 be non-unital.In an architecture with memory and readout qubits, the presence of a reset operation is essential but not on its own sufficient: as noted earlier the quantum channel must additionally contain input-dependent operations on both memory and readout qubits to prevent scrambling of memory qubits and endow the QRC state with the fading memory property.To address an important example, it is straightforward to confirm that any channel with input-dependent operations on only memory qubits U(u n ) and an arbitrary set of controlled-gates from memory to readout qubits (e.g.Fig. 5(d) in Ref. [30]) is unital on memory qubits and therefore lacks a persistent temporal memory.We note that in such cases the reservoir can still be trained to implement its function in the transient [12], but genuine online learning will not be possible.The QVT presented here prescribes how to avoid such pitfalls when designing a quantum channel for temporal data processing: one can simply check whether the resulting P 0 is a unital map.We have not carried out here an exhaustive study of the optimal design principles for more complex or general classes of tasks, but we hope that the simple and fundamental guidelines we have followed for designing an experimental reservoir to accurately carry out equalization on RF-encoded messages illustrates the utility of QVT in the design of a hardware reservoir.
By enabling online learning in the presence of losses, NISQRC paves the way to harness quantum machines for temporal data processing in far more complex applications than the CE task demonstrated here.Examples include spatiotemporal integrators, and ML tasks where spatial information is temporally encoded, such as video processing.Recent results provide evidence that the most compelling applications however lie in the domain of machine learning on weak signals originating from other, potentially complex quantum systems [15,56] for the purposes of quantum state classification.In tackling such increasingly complex tasks, the scale of quantum devices required is likely to be larger than those employed here.The NISQRC framework can be applied irrespective of device size; however, its readout features at a given time live in a K = 2 R dimensional space.For applications requiring a large R, the exponential growth of the feature-space dimension may give rise to concerns with under-sampling, as in practice the available number of shots S may not be sufficiently large.In such large-R regimes, certain linear combinations of measured features can be found, known as eigentasks, that provably maximize the SNR [9] of the functions approximated by a given physical quantum system trained with S shots.Eigentask analysis provides very effective strategies for noise mitigation.In Ref. [9] the Eigentask Learning methodology was proposed to enhance generalization in supervised learning.For the present work, such noise mitigation strategies were not needed as the size of the devices used were sufficiently small to efficiently sample.An interesting direction is the application of Eigentask analysis to NISQRC, which we leave to future work.
The present work, and the availability of an algorithm for information processing beyond the coherence time, presents new opportunities for mid-circuit measurement and control.While mid-circuit measurement is essential for quantum error correction [57], its recent availability on cloud-based quantum computers has allowed exploration of other quantum applications on near-term noisy qubits.Local operations such as measurement followed by classical control for gate teleportation have been used to generate nonlocal entanglement [58][59][60].Additionally, mid-circuit measurements have been employed to study critical phenomena such as phase transitions [61][62][63] and are predicted to allow nonlinear subroutines in quantum algorithms [64].The present work opens up a new direction in this application space, namely the design of self-adapting circuits for inference on temporal data with slowly-changing statistics.This would require dynamic programming capabilities for mid-circuit measurements, not employed in the present work.We show here that implementing even the relatively simple CE task challenges current capabilities for repeated measurements and control; having a means to deploy more complex quantum processors for temporal learning via NISQRC can push hardware advancements to more tightly integrate quantum and classical processing for efficient machine-based inference.

A. Generating features via conditional evolution and measurement
Here we detail how an input-output functional map is obtained in the NISQRC framework.The quantum system is initialized to ρMR 0 = ρM 0 ⊗ |0⟩⟨0| ⊗R , where ρM 0 is the initial state, which is usually set to be |0⟩⟨0| ⊗M .Then, for each run or 'shot' indexed by s, the process described in the following paragraph is repeated.
Before executing the n-th step, the overall state can be described as ρM,cond n−1

⊗ |0⟩⟨0|
⊗R (usually pure), where the superscript cond emphasizes that the memory subsystem state is generally conditioned on the history of all previous inputs {u m } m≤n−1 and all previous stochastic measurement outcomes.The readout subsystem state is in a specific pure state, which can be ensured by the deterministic reset operation we describe shortly.Then, the current input u n is encoded in the quantum system via the parameterized quantum channel U(u n ), generating the state ρMR,cond ⊗R .In this work, U(u n ) takes the form of continuous evolution under Eq.( 6) for a duration τ , or the discrete gate-sequence Û (u n ) depicted in Fig. 4(a).The R readout qubits are then measured per Eq. ( 2), and the observed outcome is represented as an R-bit string: b Here we consider simple 'com-putational basis' (i.e.σz ) measurements, where each bit simply denotes the observed qubit state.A given outcome j occurs with conditional probability Tr( Mj ρMR,cond n ) as given by the Born rule, and the quantum state collapses to the new state ρM,cond n ⊗ |b j ⟩⟨b j | associated with this outcome.Finally, all R readout qubits are deterministically reset to the ground state (regardless of the measurement outcome); the quantum system is therefore in state ρM,cond n

⊗ |0⟩⟨0|
⊗R .This serves as the initial state into which the next input u n+1 is encoded, and the above process is iterated until the entire input sequence u is processed.It is important to notice that ρM n depends on the observed outcome in step n − 1 and thus the quantum state and its dynamics for a specific shot is conditioned on the history of measurement outcomes {b (s) i (m)} m<n .By repeating the above process for S shots, one obtains what is effectively a histogram of measurement outcomes at each time step n as represented in Fig. 1.The output features are taken as the frequency of occurrence of each measurement outcome, as in Ref. [9]: ) counts the occurrence of outcome j at time step n.These features are stochastic unbiased estimators of the underlying quantum state probability amplitudes [9].As noted in the main text, the final NISQRC output is obtained by applying a set of time-independent linear weights to approximate the target functional ȳn = w • X(n).Importantly, during each shot s ∈ [S], we execute a circuit with depth N ; the total processing time is therefore O(N S).If instead one re-encoded N m previous inputs prior to each successive measurement the processing time is O(N m N S): N m = O(N ) if the entire past sequence is re-encoded as is conventionally done in QRC [26,27,50].

B. The quantum Volterra theory and analysis of NISQRC
At any given time step n, the conditional dependence on previous measurement outcomes, presented in Methods' subsection "Generating features via conditional evolution and measurement", is usually referred to as backaction.Defining ρMR n as the effective pre-measurement state of the quantum system at time step n of the NISQRC framework, quantum state evolution from time step n − 1 to n can be written via the maps: which describes the reset of the post-measurement readout subsystem after time step n − 1, followed by input encoding via U(u n ) into the full quantum system state.With an eye towards the construction of an I/O map, it proves useful to introduce the expansion of the relevant single-step maps U(u) and C(u) in the basis of input monomials iterative application of Eq. ( 9), ρMR n can be written as: The measured features x j (n) can then be obtained via x j (n) = Tr( Mj ρMR n ).In the Supplementary Note 3, we show that these x j (n) obtained using the NISQRC framework can indeed be expressed as a Volterra series (11) in the infinite-shot limit.The existence of this manifestly time-invariant form is only possible due to the existence of an information steady-state, guaranteed for a quantum mechanical system under measurement.
Due to fading memory, the Volterra kernel h  11)).The evolution of ρMR n upto step n − n k , namely for all i < n − n k , is thus determined entirely by the null-input superoperator P 0 .Then the existence of a Volterra series simply requires the existence of an asymptotic steady state for the memory subsystem, As shown in the Supplementary Note 3, such a fixed point is usually ensured by the map ) being a CPTP map in generic quantum systems.This immediately indicates the fundamental importance of P 0 , the operator that corresponds to the single-step map of the memory subsystem under null input: it determines the ability of the NISQRC framework to evolve the quantum system to a unique statistical steady state, guaranteeing the asymptotic time-invariance property, and hence the existence of the Volterra series.
One byproduct of computing infinite-S features {x j (n)} is that it enables us to approximately simulate { Xj (n)} in a very deep N -layer circuit for finite S, without sampling individual quantum trajectories under N repeated projective measurement described in Methods' subsection "Generating features via conditional evolution and measurement".In fact, given any n, once we evaluate a probability distribution {x j (n) ≥ 0} satisfying j x j (n) = 1, we can i.i.d.sample under this distribution vector for S shots and construct the frequency { Xj (n)} as an approximation of { Xj (n)}.The validity of this approximation is ensured by the additive nature of loss functions in the time dimension.More specifically, given )) in all orders of 1 S -expansion for any n ∈ [N ], as long as Q is large enough.This is because the probability distribution of { Xj (n)} is exactly the same as the distribution (marginal in time slice) of { Xj (n)}.Therefore, )) is a good approximation of L .
In Supplementary Note 2 and Supplementary Note 3, we show that without the reset operation, the fixed-point memory subsystem density matrix is the identity, ρMR FP = Î⊗L /2 L .While this steady state is independent of the initial state and therefore possesses a fading memory, it can be shown that the I/O map it enables is entirely independent of all past inputs as well, so that Volterra kernels h (j) k = 0 for any k ≤ 1.This yields a trivial reservoir, unable to provide any response to its inputs u.Such single-step maps C(u) are referred to as unital maps (maps that map identity to identity), and must be avoided for the NISQRC architecture to approximate any nontrivial functional.The inclusion of reset serves this purpose handily, although we have found certain improper encodings with reset to still result in unital maps C(u) (e.g., setting n T = 1 in the circuit ansatz depicted in Fig. 4(a)).
A more rigorous sufficient condition for obtaining a nontrivial functional map, referred to as fixed-point nonpreserving map in the main text, is that C(u) does not share the same fixed points for all u.It is equivalently We will prove the importance of this criteria in Supplementary Note 3. The breaking of this criteria will lead to a memoryless reservoir for all earlier input steps: if A similar result for quantum reservoirs characterized by quantum channels can also be found from Theorem 2 in Ref. [53].
C. Spectral theory of NISQRC: memory, measurement, and kernel structures Recall that we can always define the spectral problem P 0 ρM α = λ α ρM α where ρM α are eigenvectors that exist in the (2 M ) 2 = 4 M -dimensional space of memory subsystem states, and whose eigenvalues satisfy The importance of the spectrum of P 0 is obvious from the definition of ρM FP already.As ρM FP is the fixed point of the map defined by P 0 , it must equal the eigenvector ρM 1 since λ 1 = 1.Then writing the initial density matrix in terms of these eigenvectors, ρM This not only reproduces the result lim n→∞ P n 0 ρM 0 = ρM FP but also shows that the approach to the fixed point ρM FP = ρM 1 must be determined by the magnitude of λ 2 ; the smaller the magnitude, the faster terms for α ≥ 2 decay and hence the shorter the memory time.
To see more directly how the spectrum of P 0 influences memory of inputs, it is sufficient to analyze the Volterra kernels in Eq. ( 1).Focusing on single-time contributions from u n−p to x j (n) at all orders of nonlinearity (multi-time contributions are exponentially suppressed, see Supplementary Note 4), these may be expressed as which can be viewed as a spectral representation of Volterra kernel contributions to the jth measured feature obtained via POVM Mj .Here, α1 u k define 4 M − 1 internal features, so-called as they depend only on input encoding operators via αα ′ ρM α , and are in particular independent of the measurement scheme.Nontrivial F α (u) and c (k) α1 can be guaranteed if P k ρM FP ̸ = 0 for some k ≥ 1.The dependence of observables on the measurement basis is via coefficients ν Crucially, the weighting of F α (u n−p ) for p steps in the past is determined by eigenvalues λ p−1 α of P 0 .For each α ≥ 2, it vanishes when we take long time limit p → ∞.This property is usually referred as fading memory.It also clearly defines a set of distinct, but calculable, memory fading rates {|λ α |} α≥2 .
Importantly, the ability to construct Volterra kernels and internal features enable us to approximately treat the infinitedimensional function x j (n) = F j (u ≤n ) as a function with support only over a space with effective task dimension d eff = O(n M ), representing d eff time steps in the past: and we can interpret the fading memory functional as a function: In other words, at any given time NISQRC can approximate nonlinear functions that live in a domain of dimension d eff .

D. IBMQ implementation
We recall that the encoding circuit for the experimental IBMQ implementation in Results' subsection "Experimental results on quantum system" describes a composite set of single and two-qubit gates repeated n T times.
Here R x,z are composite Pauli-rotations applied qubit-wise, e.g.
i ′ } for neighboring qubits i and i ′ along a linear chain in the device and some fixed J.The rotation angles θ x,z,I are randomly drawn from a positive uniform distribution with limits [a, a + δ], where a = τ n T θ x,z,I min and δ = τ n T ∆θ x,z,I .
We find that letting the number of Trotterization steps n T = 3 is sufficient to generate a well-behaved null-input CPTP map P 0 .Our hyperparameter choices are further tuned to ensure a memory time n M commensurate with the CE task dimension.The particular hyperparameter choices for the plot in Fig. 4 are θ x,z,I min = {1.0,0.5, 0.1}, ∆θ x,z,I = θ x,z,I min , J = 1, n T = 3, and τ = 1.
In the experiment, mid-circuit measurements and qubit resets are performed as separate operations, due to the differences in control flow paths between returning a result and the following qubit manipulation [55].Related hardware complexities restrict us to a slightly shorter instance of the CE task than considered in Results' subsection "Practical machine learning using temporal data", with messages m(n) of length N = 20, submitted in batches of 200 jobs with 100 circuits each and 125 observations (shots) per circuit in order to prevent memory buffer overflows.Regardless, using crossvalidation techniques, we ensure that our observed training and testing performance is not influenced by limitations of dataset size.We also forego the initial washout period needed to reach ρ MR FP for similar reasons.Finally, the W i,i ′ (J) rotations in the two-qubit Hilbert space that implement W(J) are generated by the native echoed cross-resonance interaction of IBM backends [65], which provides higher fidelity than a digital decomposition in terms of CNOTs for Trotterized circuits [66].

DATA AVAILABILITY
The data generated for numerical results in this study have been deposited in the Github repository under accession link https://github.com/skhanCC/NISQRC-Codes[67].The raw experimental data obtained from ibmq algiers are not available in the Github repository due to its huge size, and its access can be be made available to interested parties upon request.The processed experimental data are available at the Github repository.The data of experimental parameters in this study are provided in the Supplementary Note 9.No external data was used in this study.experiments on IBM Quantum platform and analyzed the data.H.E.T supervised the project.F.H., S.A.K., G.A. and H.E.T. wrote the manuscript based on the contributions from all authors.
The underlying dynamical system we analyze in this article consists of L = M + R qubits, with M qubits serving as memory qubits and R qubits serving as readout qubits.The evolution is governed by a Hamiltonian that is linearly parameterized by a one-dimensional variable u ∈ [−1, 1] (serving as input): We choose a form of Ĥ0 and Ĥ1 that can be implemented in a quantum annealing system or analog quantum simulator in a hardware-efficient way: The coupling strength J i,i ′ , transverse x-field strength η x i = η x + ε x i and longitudinal z-drive strength η z i = η z + ε z i are pre-selected via randomness: J i,i ′ ∼ Unif[0, J max ], ε x i ∼ ε x rms × N (0, 1) and ε z i ∼ ε z rms × N (0, 1).One thing that needs to be emphasized is that the encoding scheme Supplementary Equation 1 is general enough such that encoding Eq. ( 7) is merely an illustrative example.A variety of Ĥ0 , Ĥ1 can be employed as long as they are resource-efficiently realized in a physical platform.
In theory, the domain Z of n is infinite.However in practical experiments, it is impossible to feed an input sequence from infinite past n = −∞ to infinite future n = ∞.Thus we cutoff infinity of time-step index into n ∈ The NISQRC framework then involves a continuous pipeline of evolution under a quite arbitrary superoperator U(n) (not restricted to the linearly parameterized Hamiltonian form we consider in the main text), followed by measurement, repeated until all inputs {u n } have been processed by the QRC.The inclusion of measurement with stochastic outcomes interleaved with evolution steps, as opposed to at the final step, makes our knowledge of the QRC state conditional on the entire measurement history.For example, starting from the initial state ρMR 0 and evolving under U 1 at time step n = 1, the subsequent measurement yields a measurement outcome X n = i n , where i n ∈ {0, 1} for a single readout qubit.The post-measurement state ρMR,cond 1 is then conditioned on the measurement result at time step n = 1, as indicated by the superscript cond.For an arbitrary time step n, this conditioning thus extends to the entire measurement history {X 1 , X 2 , . . ., X n−1 }.The entire pipeline can be viewed schematically as below: It is not hard to show that this process is equivalent to the quantum non-demolition scheme proposed in Ref. [2].
In practice, we are often interested not in the result of a single shot, but of the ensemble average computed over many shots; in the limit of infinite-sampling, this defines the readout features x(n) computed via ensemble averages over an infinite number of repeated shots of their stochastic conditional counterparts X n : Computing this expectation using individual measurement shots would be the standard approach in any experimental NISQRC realization, but is prohibitively expensive for this analysis.This is not least because of the dependence of X n at any time step n on the entire measurement history {X 1 , X 2 , . . ., X n−1 }, a complexity that scales very unfavourably with QRC size and the total number of time steps N .Instead, we show that the expectation can be efficiently evaluated -crucially, accounting for the conditional dynamics due to interleaved measurements -to yield a simplified expression for the infinitely-sampled readout features in terms of an effective, ensemble-averaged density matrix ρMR n , namely x(n) = Tr( M ρMR n ).To proceed, we note that, by mathematical induction, the conditional state with associated measurement record while the probability of obtaining this measurement record is simply In order to further simplify this expression, we observe the following identity for any Â ∈ C 4×4 , which can be verified by direct computation i=0,1 where the matrices Î = 1 0 0 1 and 1 = 1 1 1 1 , and the notation • represents the Hadamard product (element-wise product): Supplementary Equation 9 enables us to introduce the measurement-induced decoherence superoperator M: Therefore, according to Supplementary Equation 8, the unconditional expectation of the random variable X n can be computed by contraction: where we used Supplementary Equation 10 and M = in i n Pin .This expression naturally leads to the identification of the term in square brackets as the effective density matrix at time step n, ρMR , such that computing the trace with respect to this density matrix provides any readout feature at time step n in the infinite sampling limit, x(n) = Tr( M ρMR n ).The generalization to a QRC with L = M + R ≥ 3 and input sequence {u n } is now straightforward: U n is replaced with U(u n ), while the measurement-induced decoherence superoperator M generalizes to: With these changes, the effective density matrix at time step n for the NISQRC framework without reset is given by Note that ρMR n accounts for both any time-dependent unitary dynamics via U n , as well as the role of repeated measurements via recurrent applications of M.

a. Thermalization induced by repeated measurements without reset
We need to point out that even if the circuits have similar structures to those used in measurement-induced phase transition [3]: at step n associated with unitary evolution U n , qubits indexed by a random subset I n ⊆ [L] will be measured.In this scenario, the effective state evolution is similar ρMR ρMR 0 , the only difference is that measurement-induced decoherence superoperator M n now is no longer a time-independent map For any overall state ρMR , the Frobenius distance ρMR − I ⊗L 2 L 2 F will never increase after either unitary evolution U or measurement M I : where the proof employs that fully mixed state Î⊗L 2 L is the simultaneous fixed point of U and M I (equivalently, both maps are unital CPTP map).The non-increasing purity implies that The final QRC state therefore has no memory of the initial state ρMR 0 .As a result, in previous works [2,4] this type of evolution has been employed to equip QRCs with the fading memory property.However, note that the final state is also entirely independent of the input u(n), which renders it incapable of performing any useful computations on this input.Hence input-dependent unitary evolution combined with readout only does not yield a useful QRC.We show next how a simple modification of the measurement protocol can allow fading memory without yielding a trivial I/O map.

Quantum dynamics under measurement and reset
For notational simplicity, we once again analyze a system with M = 1 memory qubit and R = 1 readout qubit (namely L = M + R = 2).We apply Pauli z measurement on the readout qubit at each QRC step, the corresponding observable is M = Î ⊗ |1⟩⟨1|.Since now we apply the conditional reset.The measurement process is described by a POVM measurement (i = 0, 1): Also, for more general M > 1 and R > 1 we used in the main text, we can still introduce the effective density matrices ρMR n in NISQRC having the same expression where C(u)ρ M = Tr R U(u) ρM 0 ⊗ |0⟩⟨0|

⊗R
. Hence, we finish deriving the expression of ρMR n .

Tr(P
If we take P k ρM = Therefore, we conclude a very useful property that for any memory density matrix ρM and any k ≥ 1.
2. R k and P k for linear Hamiltonian encoding scheme by regrouping the BCH formula We now evaluate the u-expansion of U(u)ρ MR = e −iτ Ĥ(u) ρMR e iτ Ĥ(u) .Central to this expansion is the Baker-Campbell-Hausdorff (BCH) formula, which allows us to write this expression in the series form Using the explicit form Ĥ(u) = Ĥ0 + u Ĥ1 , we can compute the superoperator coefficient of any term in the series: . . .Note that each term in the series can be viewed as a series in u k instead.Furthermore, each appearance of u Ĥ1 in Ĥ(u) contributes exactly one factor of u.This allows us to determine the coefficient of u k in the qth term: Here, the summation is over all q k possible combinations { Ĉ1 , Ĉ2 , • • • , Ĉq } which is an ordered set with k instances of Ĥ0 and (q − k) instances of Ĥ1 .This expression allows us to regroup the BCH formula not by the parameter q as in Supplementary Equation 34, but by powers u k of the input.We therefore arrive at the desired form of Supplementary Equation 26, with 3. Functional I/O map: time-invariance and Volterra kernels Our work in the previous subsection allows us to express the action of individual superoperators U(u) and C(u) on a general ρMR as a u-expansion at every time step.The dynamical map defined by our time-dependent NISQRC framework involves the repeated application of these superoperators for distinct inputs u n , so that the output at time step n may have a complicated dependence on prior inputs u ≤n .We are now in a position to extract this dependence explicitly.To do so, we simply substitute our u-expansions for the superoperators U(u) and C(u) into the evolution equation Supplementary Equation 25defining ρMR n at an arbitrary time step n, i.e.
Then, the density matrices at time step n attain the formal expression: Before evaluating the readout features x j (n), we need to simplify Supplementary Equation 38as much as possible.The starting point is first looking at the simplest contribution from term u n−1 to ρMR n (namely the one-step backwards linear contribution).This means that we can let k The associated prefactor is Similarly, analyzing contribution from term u n to ρMR n+1 (that is, let it ensures the approximation and hence Supplementary Equation 39 and Supplementary Equation 40 are asymptotically the same.Such property is usually referred as (asymptotic) time-invariance.In fact, we can further weaken this requirement that all peripheral spectrum λ α (namely those eigenvalue with magnitude |λ α | = 1) are λ α = 1.For example, for a fully connected quantum reservoir with M +R qubits, if J i,i ′ are constant for every coupling pair and η x i , η z i are also constant for every qubit, then the numerical results show that the fixed points of ρM 0 will have a degeneracy of Catalan numbers (2M )! M !(M +1)! .In this case, the fixed point lim n→∞ P n 0 ρM 0 = ρM FP still exists but will depend on initial state ρM 0 .The above calculation works for any contribution terms in ρMR n .This establishes all analytical expressions of Volterra series kernels.The leading order kernels can be written down compactly: • The zero-th order Volterra kernel: • The first order Volterra kernel (n 1 ≥ 0): illustrative for this Supplementary Note to mostly be concerned with a single past time step's contribution.To be more specific, if we focus on the contribution from u n−p to x j (n) (where p ≥ 1).For this history record contribution, Especially, if α = 1, then c (k) α1 = 0 for any k ≥ 1, according to Supplementary Equation 32.That is why the summation over α starts from α = 2, and it only gives us 4 M − 1 internal features (see Supplementary Figure 4 as an example).

Supplementary Note 5: Relation between functional-independence and Jacobian rank
In this Supplementary Note we analyze the functional-independence of readout features in the NISQRC framework.Assuming a finite-dimensional input space u = (u 1 , u 2 , • • • , u N ), then NISQRC readout features define K finite-dimensional functions (assuming K ≤ N ), x k (u 1 , u 2 , • • • , u N ), k ∈ {0, • • • , K − 1}.An important question is whether these K functions are in fact functionally-independent from one another, since their inter-dependence can impose a limitation on their usefulness for functional approximation using the NISQRC framework.
If the K functions are functionally-dependent, namely there exists some K-variate function G such that: for all non-negative integer p ∈ N.For generic input sequence {u −∞ , • • • , u n−1 , u n }, there doesn't exists such observable K−1 j=1 c j Mj such that expectations of ∂ ρMR n ∂un−p for any p ∈ N always vanish, which is a contradiction.This results shows that in principle, the linear combination of quantum probability readout will yield a function family whose gradient space is much more abundant, because usually the feature number K − 1 is much larger than the readout qubit number R.

k
(n 1 , • • • , n k ) characterize the dependence of the systems' measured output features at time n on its past inputs u n−nκ .Hence the support of h (j)

Fig. 2 (
b) for various values of τ .The unique eigenvector corresponding to the largest eigenvalue λ 1 = 1 is special, being the fixed-point of the memory subsystem, ρM 1 = ρM FP , reached once transients have died out.

FIG. 3 .
FIG. 3. Numerical results for the channel equalization (CE) task with Hamiltonian ansatz.(a) Error rates on test messages for the CE task with a Hamiltonian ansatz (2 + 4)-qubit QRC for two distinct connectivities shown in (b) The fully-connected QRC in red has Jacobian rank RJ = 2 R − 1 = 15 and is shown for both S → ∞ (circles) and finite S = 10 5 (⋆), whereas the split QRC has RJ = 2(2 2 − 1) = 6 and only S → ∞ is plotted in magenta.These are compared with the error rates of naive rounding (black dash-dots) and logistic regression on the current signal (yellow +, see Supplementary Note 6), and the exact channel inverse (blue dashed).(c) Performance of connected QRC on SNR = 20 dB test signals (solid) of increasing length Nts ≤ 5000, with shots S = 10 5 .Training error on N = 100-length messages is indicated for comparison in dashed lines.Without reset (red) or using 4 ancilla qubit ansatz with quantum non-demolition (QND) readout (proposed in Ref.[30], green), the algorithms both fail, approaching the random guessing error rate and showing that both architectures suffer from the thermalization problem.Performance is only slightly reduced from the dissipationfree case (blue) when strong decay T1 = 10τ is included (purple).All error rates in (c) are averaged over 8 different test messages.

22 FIG. 4 .
FIG. 4.Experimental results for the channel equalization task with circuit ansatz.(a) (3 + 4)-qubit linear chain of the ibm algiers device used to perform the CE task.Filled colors represent qubit T1 time according to the displayed colorbar, for the specific experimental run with the split chain.Qubits indexed {8, 14, 19} are used for memory and qubits {5, 11, 16, 22} for readout, and gatedecomposition of the encoding unitary Û (un) is depicted.Removing gates shaded in brown yields two smaller chains to explore the role of connectivity, while removing reset operations (shaded peach) allows switching from a non-unital to a unital I/O map.(b) Testing error rates for the SNR = 20 dB CE task of Results' subsection "Practical machine learning using temporal data" with N = 20 on the ibm algiers device in filled circles and in simulation in open circles, as a function of number of shots S. The connected circuit in blue outperforms the split circuit in brown and the circuit without reset in peach.For comparison, we plot the testing error rate of logistic regression (yellow line), as well as random guessing (black dashed line).
(j) k (n 1 , • • • , n k ) characterizes the dependence of the systems' output at time n on inputs at most n k steps in the past (recall n 1 ≤ • • • ≤ n k , see Eq. ( 0 and k n = 1) gives associated prefactor R 0 P 1 P n−1 term Supplementary Equation 39 and Supplementary Equation 40 are analytically different.However, with the existence of fixed point state lim n→∞ observable Mj 's response to different internal features F α (u) where F α (u) = u) ∂x1 ∂u N (u) • • • ∂x K−1 ∂u N (u) 0 , x 1 , • • • , x K−1 ) ∂G ∂x1 (x 0 , x 1 , • • • , x K−1 ) . . .∂G ∂x K−1 (x 0 , x 1 , • • • , x K−1 ) gradients ∇ u x 0 (u), ∇ u x 1 (u), • • • , ∇ u x K−1(u) must be linearly dependent at all points.Therefore, if {x j (u)} j∈[K] are functionally-dependent, then the gradients {∇ u x j (u)} j∈[K] must be linearly-dependent.Equivalently, it suffices to prove the functional-independence of {x j (u)} j∈[K] by showing that {∇ u x j (u)} j∈[K] are linearly-independent at almost all points u.Now we argue by contradiction that K − 1 gradients of readout features x j (n) = Tr Mj ρMR n are functionally-independent if there is no particular symmetry in the reservoir.We first select { Mj } as the moment representation to remove the trivial functional dependence that their summation is constant.Suppose there exists coefficients c 1 , c 2 , • • • , c K−1 such that K−1 j=1 c j ∇F j (u ≤n ) = 0. Notice that ∂x j ∂u n−p c j ∇F j (u ≤n ) = 0 implies that is the most different one since its associated matrix ρM 1 = ρM FP will remain unchanged when applied by P 0 while other modes decay to zero: lim n→∞ P n 0

Table 2 .
Device parameters for split QRC with mid-circuit measurement and deterministic reset (brown line in Fig. 4(b) of the main text).